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ABSTRACT 

We use the Cambridge stellar evolution code STARS to model the evolution 
and nucleosynthesis of zero-metallicity intermediate-mass stars. We investi- 
gate the effect of duplicity on the nucleosynthesis output of these systems and 
the potential abundances of the secondaries. The surfaces of zero-metallicity 
stars are enriched in CNO elements after second dredge up. During binary 
interaction, such as Roche lobe overflow or wind accretion, metals can be re- 
leased from these stars and the secondaries enriched in CNO isotopes. We in- 
vestigate the formation of the two most metal poor stars known, HE 0107-5240 
and HE 1327-2326. The observed carbon and nitrogen abundances of HE 0107- 
5240 can be reproduced by accretion of material from the companion-enhanced 
wind of a 7 Mq star after second dredge-up, though oxygen and sodium are 
underproduced. We speculate that HE 1327-2326, which is richer in nitrogen 
and strontium, may similarly be formed by wind accretion in a later AGB 
phase after third dredge-up. 

Key words: binaries: general, stars: abundances, stars: AGB and post- AGB, 
stars: carbon, stars: evolution, stars: individual: HE 0107-5240, HE 1327-2326 



1 INTRODUCTION 

Mass loss from primordial stars is very important because any metals released from them 
enrich the interstellar medium(ISM), increasing the metallicity of subsequent generations 
of stars. However mass loss from single stars, though not well understood, is thought to 
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be small compared to solar metallicity stars. To be efficient radiation-driven stellar winds 
require metal lines (IMarigo. Chiosi fc Kudritzkill2003l ). Thus the winds of zero-and very low- 
metallicity stars may be very weak, owing to the absence of metals at the surface. However 
we do not know how the mass loss depends on metallicity during the asymptotic giant 
branch (AGB) phase, so it is very uncertain how much this is affected by the low metallicity. 
On the other hand, mass loss from stars during binary interaction is more well-defined and 
understood. If the mass transfer is not conservative, the surroundings are enriched in metals. 
If some of the material lost by the primary is deposited on to the secondary then its surface 
abundance becomes CNO rich but otherwise metal-poor. This could be one way to form 
the observed carbon-enhanced metal-poor stars and, at the same time, make a significant 
contribution to the contamination of the ISM. Observed carbon-enhanced metal-poor stars 
have a large ove rabundance of carbon, [C /Fe] > +1.0, and they constitute 20% of stars with 
[Fe/H] ^ —2.0 (jBeers &: Christliebl 120051 ) . This suggests a significant amount of carbon was 
produced in the early Universe. Binary interaction is one possible way to create carbon-rich 
stars. Alternatively, it is also p ossible that some carbo n-rich stars formed from material 
enriched by type-II supernovae (jUmeda fc Nomotoll2003l ). 



1.1 Possible formation of binary intermediate-mass primordial stars 



The primordial, first generation of stars, commonly referred to as Population-Ill stars or zero- 
metallicity stars, should have the composition of the ISM just after Big Bang nucleosynthesis 
and hence should be devoid of metals. It has been a popular belief that, in the absence of 
heavy elements and dust grains, co oling mechanisms are inefficient and favour the formation 



of massive or very massive stars ( IBromm. Copp: 



fc Larson 



3) 



1999) because the Jeans mass 



would be large. However it has also been shown ( IPalla. Salpeter fc Stahlerlll983l ) that even 
a small fraction of molecular hydrogen can provide a significant contribution to the cooling 
via rotational and vibrational transitions. The resulting Jeans mass of a pure H and He 
cloud is then relatively small and may even fall below 0.1 M & . Owing to this complexity 



and our lack of understanding of star formation, the initi a 



metallicity stars remains uncertain. 



Nakamura fc Umemura 



mass function (IMF) of zero- 



(j200lh suggested a bimodal IMF 



with peaks around 1 M & and 100 M . Hence, intermediate-mass zero- metallicity stars should 
be investigated. 
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Saigo. Matsumoto fc Umemural (12004 ) explored the formation of Population-Ill binaries, 



by simulating the collapse of a rotating cylinder with a three-dimensional, high-resolution 
nested grid. They found that low angular momentum is key to the formation of binaries. If 
the initial angular momentum is small enough the runaway collapse of the primordial gas 
cloud stops at some point and forms a rotationally supported disc instead of continuing to 
collapse into a single core. After accreting the surrounding cloud, the disc can eventually 
undergo a ring instability and fragment into a binary. However, if the initial angular momen- 
tum is relatively large, a bar-type instability arises and the cloud collapses to a single star 
through rapid angular momentum transfer. Their results showed that a significant fraction 
of Population-Ill stars is expected to form in binary systems, regardless of their mass. Hence 
the effect of binaries on the evolution and nucleosynthesis of primordial stars can be very 
important. 



2 THE STARS CODE 



We use the Cambridge stellar evolution code STARS to calculate t he evoluti o n of p rimordial 



Eggletonl (119711) and ha s 



Pols et al. 



dl995h . 



intermediate-mass stars. This program was originally written by 
been updated by many authors. The last major update was undertaken by 
It solves the equations of stellar structure, nuclear burning and mixing simultaneously, unlike 
most codes which treat mixing in a separate step. Another unique fe ature of the STARS code 



is the use of a self-adaptive non-Langrangian, non-Eulerian mesh. 



Eldridge fc Tout 
J 



updated the opacity tables to use the latest OPAL opacities ( llglesias fc Roger i H996I) and 
these new tables also account for the change in opacity with variations in the carbon and 
oxygen abunda nces. We use their zer o-metallicity opacities in this work (Eldridge, private 



communication; 



were recently resurrected by 



Ferguson et al.l 120051). The nucleosy nthesis subroutines of the STARS code 



Stancliffe et al. 



(120051 ). These make use of the output from 



the structure code to determine the mixing and nucleosynthesis. In addition to the seven 
main isotopes that are followed in the structure code ( 1 H, 3 He, 4 He, 12 C, 14 N, 16 O, 20 Ne) the 
evolution of 37 other isotopes, including all stable and a few unstable isotopes from 2 D to 
34 S and the major iron group elements, are currently modelled by the nucleosynthesis code. 
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3 EVOLUTION OF SINGLE STARS 



We have modelled 3 — 7 M & stars from the pre-main sequence to the beginning of thermally- 
pulsing AG B using 999 meshpoi nts without any convective overshooting. The mixing-length 



parameter (jBohm-Vitense 



19581 ) a is 1.925 based on calibration to a solar model. The he- 



lium mass fraction is chosen to be 0.25 to reflect the prediction of primordial helium abun- 
dance from the observed deut erium abundance, baryon de nsity and a spectroscopic sample 



of extragalactic H-II regions (jFukugita fc Kawasaki! 120061 ). The e volution of our mod els is 

(1200 If ) and 



in qualitative agreement with the previously published work of 



Chieffi et al. 



Siess. Livio &: Lattanziol (120021 ). even though they used a different value of the initial helium 



abundance. 

We now briefly describe the overall characteristics of the evolution of intermediate mass 
zero-met allicity stars. The evolution of zero-metallicity stars differs significantly from that 
of population-I and -II stars because of the absence of carbon, nitrogen and oxygen. Hy- 
drogen cannot burn through the CNO cycle so it burns only by the proton-proton chain. 
This is less temperature dependent so zero-metallicity stars are considerably hotter in the 
nuclear-burning regions and throughout most of the star. Their main-sequence lifetimes are 
much shorter and the burning is much less centralised. Furthermore the temperature at the 
core becomes hot enough that carbon is slowly produced from helium through the triple-a 
reaction. At this time there are protons at high temperature so the products of helium burn- 
ing end up mainly as nitrogen in CNO equilibrium. When the central carbon mass fraction 
reaches about 10~ 12 hydrogen burning switches to the CNO cycle. For stars more massive 
than 1.8 M & , the core is hot enough for helium to burn non-degenerately when central hydro- 
gen is exhausted. This mass, which is the dividing line between low and intermediate-mass 
stars, is much lower than for solar-metallicity stars because zero-metallicity stars are much 
hotter and hence have less degenerate cores. The intermediate-mass stars remain on the blue 
side of the Hertzsprung-Russell (H-R) diagram throughout central helium burning and so 
the red-giant phase and first dredge up are absent. The surface abundances are unchanged 
during the core He-burning phase. The upper mass limit for AGB stars at Z = is between 
7 M & to 8 Mq because our 8 M & model starts carbon ignition at the core before thermal 
pulses. 

When helium is depleted in the core the stars undergo a large expansion and move 
rapidly to the red side of the H-R diagram. Hydrogen and helium burning are active in two 
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1.59 x 10" 2 


1.14 


(5.0 


0.382 


9.38 x 10 
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-8 
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5.42 
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0.382 
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-9 
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-8 


3.85 x 10- 
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1.16 x 10 2 
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Table 1. Surface abundances and abundance ratios by mass fraction at the end of the early AGB before thermal pulses. 



separate shells. The energy generated by the helium shell causes expansion and cooling of the 
outermost layers and the convective envelope deepens. This is second dredge-up, although it 
is actually the first episode of dredge-up that occurs in zero-metallicity stars above 1.8 M Q . 

Enough carbon is produced that hydrogen is mainly burned through the CNO cycle in 
the H-burning shell and nitrogen is produced therein. Abundance profiles for a 7M star 
are shown in Figure HJ When the convective envelope deepens both carbon and nitrogen are 
dredged to the surface. Lower-mass stars dredge up less of the CNO isotopes because their 
convective envelopes do not penetrate so deeply and the CNO abundance in the H-burning 
region is lower because the star is cooler. The surface C/N ratio (shown in Table 1) is in CNO 
equilibrium at the temperature of the burning shell. Around 1/100 for lower- mass stars, it 
increases rapidly from 5M to 7M . This is because second dredge- up reaches down to 
the helium-burning shell for higher- mass stars. In addition, the ratio of 12 C/ 13 C increases 
from 3.5 to 1190. This requires an additional source of 12 C that is not present in the CNO 
H-burning shell but is also explained by this deeper dredge up. Because of the dredging up 
of materials from the He-burning shell, the carbon abundance increases from 4.07 x 10~ n 
to 2.64 x 10~ 6 for stars from 5M to 7M . Oxygen follows a similar trend, increasing from 
3.56 x 10~ n to 3.85 x 10 -9 . Nitrogen is already significantly enhanced at 2.56 x 10 -9 for a 
5M Q star and only increases a little more to 2.28 x 10~ 8 in the 7M model. 



4 MASS LOSS FROM AGB STARS 



It is generally believed that ma ss loss from zero-me tallity stars is very slow because radiation 
driven winds are much weaker (IMarigo et al.ll2003l ). However, the mass-loss mechanism dur- 
ing AGB evolution is not well understood. It has been suggested that the increase of the mass- 
loss rate during the AGB is caused by the onset of radial p ulsations ( jVassiliadis fc Wood 



1993 



Bloecker 



1995T ) or magnetic fields (IBlackman et al.ll200ll ) . Whether the mass loss would 



then be dependent on metallicity is not known. In addition the lower opacity of a metal- 
deficient atmosphere generally leads to more compact stars so the mass-loss rate is expected 
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Figure 1. Abundances by mass fraction profile of the inner part of a 7Mq model near the end of core helium burning. 



to be lower in any case. However, after dredge up, the surfaces o f the stars are also enriched i n 
carbon and rates may increase again. Reimers' mass-loss rate (IKudritzki fc Reimerslll978l ). 



M = 4.0 x lO- n V {L/ ^ R J R&) M & yr-\ (1) 
' (M/M & ) w 

where L is the luminosity of the star, M its mass and R its radius, repro duces the early 



AGB mass loss of higher metallicity stars reasonably well when 0.4 < rj < 3 fjStraniero et al 



20031 ). We use r\ of both 1 and 4 (which almost certainly overestimates the mass loss from 



Z = stars) for 3M , 5M and 7 M & models to simulate the mass loss during the early 
AGB phase. We define the yield of isotope i by 

Yt= [ T M(t)Xi(t)dt, (2) 
J o 

where M(t) is the mass-loss rate, T is the lifetime of the star and Xi(t) is the surface 
abundance at time t of isotope i. 

As can be seen in Tables 2 and 3, the yield is directly proportional to r\ because changing 
the value of rj has little effect on the stars' evolution because their mass-loss rates remain 
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Initial Mass/M0 


Mass Lost/ik/Q 


Yl2 C /M Q 


Yu n /Mq 


Yib /Mq 


3.0 


0.077 


1.71 x 10~ 17 


1.94 x 10" 15 


5.05 x 10~ 17 


5.0 


0.052 


1.22 x 10~ 12 


6.95 x 10" 11 


9.07 x 10 -13 


7.0 


0.034 


4.82 x 10" 8 


3.61 x 10" 10 


6.77 x lO" 11 


Table 2. Mass loss and yield of CNO before thermal pulses with r] = 1 for 


single star. 




Initial Mass/M Q 


Mass Lost/M 


Yi2 C /M Q 


yi4 N /M Q 


Yie /M @ 


3.0 


0.32 


8.01 x 10~ 17 


9.06 x 10 -15 


2.36 x 10~ 16 


5.0 


0.22 


4.79 x 10" 12 


2.74 x 10" 10 


3.63 x 10" 12 


7.0 


0.14 


1.94 x 10" 7 


1.47 x 10~ 9 


2.70 x lO" 10 


Table 3. Mass loss and yield of CNO before thermal pulses with r\ = 4 for 


single star. 





low. Our 3 M Q star releases a negligible quantity (of the order of 10~ 15 ) of metals because the 
dredge-up is not deep enough to bring many to the surface. Our 7M star has significantly 
higher yields of CNO elements, particularly carbon (4.82 x 1O~ 8 M for n = 1) and to a 
lesser extent oxygen (6.77 x 10~ n M ). This is because the convective envelope reaches very 
close to the helium burning shell to dredge up carbon and this is only possible if the star 
is more massive than about 5M . The nitrogen yield for a 5M star is 6.95 x 10 -11 M 
compared with 3.61 x 1O~ 1O M for 7 M & . A little carbon is also present outside the helium- 
burning shell because it was produced at the end of main sequence by the triple-a reaction 
when the core was hot and extended. Most of the carbon produced then was immediately 
converted to nitrogen by the CNO cycle. So in order to dredge up nitrogen the convective 
envelope needs only to reach the region in which the CN cycle was active. This first happens 
at around 5 M . Given the uncertainty in the initial mass function for zero-metallicity stars, 
we cannot make an accurate prediction of the relative contribution of different mass ranges. 
However, because the contribution of carbon from a 7M star is more than 1,000 times 
that of a 5 M star, we can see that the highe r-mass AGB st ars can dominate the carbon 



r 



contribution for a normal IMF such as that of ISalpeterl (119551 ). Similarly, the oxgyen yield 
is dominated by the 7M stars because it is nearly 100 times greater than that of a 5M 
star. Unless the ratio of 5M to 7M stars is a few times higher than Salpeter's ratio the 
nitrogen contribution from 5 M stars is negligible compared to that of 7M stars. 



5 BINARY STARS 

We expect mass transfer by Roche lobe overflow in low-metallicity binary stars to proceed 
in a similar way to higher metallicity systems. If the mass transfer is non-conservative the 
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surroundings are also enriched in metals but the degree to which mass transfer is non- 
conservative is not known. If the system is a very close binary, mass transfer begins on the 
main sequence. Otherwise, because the radius of a zero-metallicity intermediate-mass star 
only increases substantially after it leaves the main sequence, it is more likely for the mass 
transfer to begin in the early AGB phase than before. On the main sequence the star's radius 
increases by a factor of two while on the early AGB it increases by a factor of one hundred. In 
this case, when the primary fills its Roche lobe, any mass lost contains a significant amount 
of metals because dredge up has already occured. In cases when mass transfer begins before 
second dredge up, such as the Hertzsprung Gap, CNO isotopes can still be transferred 
because the CNO-rich H-burning shell is eventually exposed at the surface. 

We have made several binary models with primary masses ranging from 5 to 7 M & and an 
initial mass ratio of two-thirds, with different periods (see table 4). These models are evolved 
without any tidal interaction. The stellar evolution of the secondary is also not modelled 
because a lower mass secondary would not have evolved when mass transfer starts. The 
secondary star is thus treated as a point mass. The variation in the period of the system is 
modelled by considering the change in angular momentum during binary interaction. Mass 
transfer begins when the primary is in the Hertzsprung Gap for systems shown in the table. 
Longer periods (not listed) could result in a common envelope because mass transfer starts 
when the primary's envelope is deeply convective. If tides were included, the periods would 
decrease because orbital angular momentum would go to spin up the expanding primary. 
We do not explicitly model tidal effects so our models really correspond to slightly larger 
initial periods. In the case of a 7M primary star the envelope can be completely lost to 
leave a 1 M Q CO white dwarf. Even when 6 M & is transferred to the companion, we are 
only interested in the mass released from the primary, and so treat the companion as a 
unevolved point mass. This does not have any significant effect on the result. Provided that 
the mass transfer rate remains low, the radius increase due to accretion by the secondary 
would not severely disrupt the binary system. Significant masses of metals are transferred 
from the primary, up to 4.00 x 1O" 5 M , 5.23 x 1O~ 7 M and 1.93 x 1O" 7 M of C, N and 
O respectively. Further increases in the orbital period would only have a small effect on the 
compositions of the material lost because mass transfer starts after the metals have already 
been dredged up. As the binary period increases the remnant mass slightly increases but 
this has only a very small effect on the total amount transferred from the primary. 

The mass of metals transferred from the primary star is more than one hundred times 
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higher than the mass of metals released during typical single-star mass loss at the same 
evolutionary point. In particular the carbon transferred from the primary envelope is nearly 
1,000 times greater. This is because the whole envelope of the primary is lost during binary 
interaction while a very small fraction is lost from a single star by the AGB wind. We have 
computed a single 7 M Q star model through all its 600 thermal pulses or so and found that 
these are too weak to lead to any third dredge up. Carbon is converted to nitrogen because of 
hot bottom burning. Later pulses are weaker and they cease altogether when the core mass 
reaches 1.1 M & . The stars then enter a non-eventful evolutionary phase while the hydrogen 
and helium burning shells grow outward. Eventually, carbon ignites degenerately at the 
centre when the core mass is 1.36 M . The subsequent thermonuclear runaway (similar to a 
type la supernova) would release sufficient energy to blow the whole star apart. At this point, 
the surface carbon and nitrogen abundances are 5.5 x 10 -7 and 2.7 x 10" 6 . The composition 
of the lost material would be very different to the binary case discussed previously, with up 
to 1.6 x 1O" 5 M of nitro gen being ejected while the yield of carbon would be of the order of 
3.3 x 10~ 6 M Q depending on mass-loss rate. This is over thirty times more nitrogen. Despite 
the envelope, we would expect the explosion mechanism of this star to be very similar to a 
Type la supernova, so the carbon and oxygen ejected from the supernova later on exceeds 
the yield from the envelope (Lau, Stancliffe & Tout, in prep). 

We also model a 7 M & primary star with a 4 d period with different mass transfer efficien- 
cies, 0.25, 0.5 and 0.75, where the mass-transfer efficiency is the ratio of the mass transferred 
to the seondary to the mass lost from the primary. A mass transfer efficiency of 0.75 means 
three-quarters of the mass lost from the primary is accreted by the secondary and the re- 
maining quarter is lost from the system. We calculate the masses of metals released from 
the binary system and list them in table 5. The increase in the metals released is already 
significant even when the mass-transfer efficiency is 0.75. The carbon, nitrogen and oxygen 
yields are 100 to 1,000 times more than the single star yield of a similar age. 

For the case of mass-transfer efficiency of 0.25, mass transfer proceeds very rapidly be- 
cause of the greater mass ratio throughout owing to the secondary accreting less material. 
Because mass is transferred from the more massive primary to the less massive secondary, 
the Roche lobe radius shrinks, while the radius of the primary increases as mass is lost from 
its convective envelope. This leads to positive feedback and further increases the mass-loss 
rate. In practice, this probably leads to common envelope evolution. The common envelope 
surrounding both stars is likely to be lost from the system, so we assume the remaining 
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Primary Mass/M Period/d Final Mass/M Y li c / M Q Y L^/ M Q Y le / M @ 
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Table 4. The masses of metals transferred from the primary Y'in various binary interactions. Some is accreted by the companion 
and rest lost to the ISM 
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Table 5. The masses of metals released from a binary system with initial primary mass of 7,Mq, q = 2/3 and period of 4 
days. 

primary's envelope is released from the binary soon after the common envelope forms. In 
this case, the total amount of CNO released is even higher because most of the primary's 
metal-rich envelope is released from the system. In table 5, we assume all the envelope of 
the primary is lost from the system at the onset of the common envelope evolution. 



6 CARBON ENHANCED METAL POOR STARS 



There have been two di scoveries of carbon-enh anced hyper-metal-poor sta rs (—6.0 ^ |Fe/H] ^ 



5.0), HE 0107- 5240 (IChristlieb et al. 



2002 



Aoki et al 



2004) and HE 1317-2326 (IFrebel et al. 



2005; 



20061 ). They are the lowest-metallicity stars observed so far. The important ele- 
mental and isotopic abundances are summarized in Table 6. 

HE 0107-5240 is not significantly enriched in strontium while the other has an unexpect- 
edly high strontium abundance. Carbon is very much enhanced relative to the solar C/Fe 
ratio, by a factor of several thousand. Nitrogen and oxygen are also strongly enhanced, being 
a factor of more than 100 above the solar N/Fe and O/Fe ratios. The C/N abundance ratio 



Abundances 


HE 0107-5240 


HE 1317-2326 


7 Mq early AGB models 


[Fc/H] 


-5.3 


-5.45 


-5.26 


[C/Fe] 


3.7 


3.90 - 4.26 


2.73 


C/N 


40 - 150 


3 


239 


[O/Fe] 


2.3 


< 4.0 


0.1 


12 C/ 13 C 


> 50 


> 5 


14200 


[Sr/Fe] 


< -0.52 


0.91 - 1.24 


Sr not enhanced 



Table 6. ID atmospheric analysis abundances and surface abundances of observed CeHMP stars and our 7Mq star with 
Z = 10 -7 after second dredge-up. 
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of 40 to 150 for HE 0107-5240 is similar to the surface abundances of the early AGB phases 
of our 6 — 7 M Q stars with Z = 10~ 7 , with the initial isotopic ratios taken as solar 

Because there is already iron enrichment, the progenitors of these stars have not evolved 
directly from Z = stars. We therefore carried out a new set of calculations at Z = 10 -7 . For 
these intermediate mass stars this does not introduce enough CNO isotopes to qualitatively 
alter their evolution from those of Z = but we model a new set at Z = 10~ 7 for consistency. 
The observed 12 C/ 13 C ratio, which is greater than 50, is not inconsistent with the models. 
We do not include strontium in our code but we can deduce that strontium is not enhanced 
in the early 7 M & AGB stars because of the absence of a neutron source. One possible 
formation scenario for HE 0107-5240 is that it is the secondary star of a wide binary system, 
of which the original primary star has already evolved into an unsee n white dwarf. This i s 
very similar to the best explanation for Ba stars of solar metallicity (IKarakas et al.ll2000l ). 
Mass was tran sferred to the seco ndary through wind accretion during the early AGB phase 

(120041 ) proposed an alternative scenerio in which the primary star 



of the primary. 



Suda et al. 



has an initial mass between 1.2 and 3M . The major difference is that, in their lower mass 
models, the controversial helium-flash-driven convective mixing is responsible for bringing 
metals to the surface, while the well-established second dredge up alone is enough to increase 
surface carbon and nitrogen in our 7 M & models. 

In order to explain the non-detection of radial velocity variation, the binary has to 
be wide enough to avoid Roche lobe overflow and common envelope evolution. Otherwise, 
the final period would be low enough for the companion to be detectable. The minimum 
initial period of the binary system must have been at least 600 d. Because of the loss of 
angular momentum in the wind, such a period would increase to at least 30 yr at the present 
day. Moreover in our model the current radial velocity of the secondary would be about 



6.5 km s . This theoretical value with the 30 



variation in radial velocity of HE 0107-5240 (jBessell. Christlieb fc Gustafsson 



yr period can explain the non-detect ion of 



3). 



2004). The 



maximum velocity variation would be 4.8 km s -1 after 4 years of observation with radial 
velocity variation of up to 13kms _1 in 15 yr. Considering the orientation of the orbit, we 
might expect a velocity variation of about lkms" 1 after 4 years and 5kms _1 in 15 yr. 



For systems that are too wide for Roche lobe overflow, iTout fc Eggletonl (119881 ) suggested 
that mass loss may be tidally enhanced by the presence of a moderately close companion. 
Their enhanced mass loss formula is 
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M = M R {1 + B max 



1 R^ 6 



(3) 



where Mr is the original Reimers mass-loss rate for a single star, R^ is the Roche-lobe 
radius and B ~ 10 4 . A fraction of the matt er may be accreted by the secondary through 



Bondi-Hoyle accretion (IBondi fc Hoyle 



194^). The system widens owing to loss of angular 



momentum and can remain too wide for Roche lobe overflow. In this case, mass transfer 
is not conservative and metals are released from the system. Because wind accretion only 
increases the secondary mass by a small amount and mass loss via winds increases the bi- 
nary separation, this mechanism provides a possible formation scenerio for carbon-enhanced 
metal-poor stars with masses around 1 M Q and no observed binary companion. These may 
be long period binaries. Wind accretion can be significant if the binary is close enough that 
the primary radius is close to the Roche Lobe radius during the e arly AGB phase . 



Using the enhanced mass- loss rate and Bondi-Hoyle accretion (IBondi fc Hoyle! 1 19441 ). we 



estimate the amount of material that can be accreted from a 7M m star. The mass accreted 



by the secondary (IHurley. Tout fc Pola 120021 ) 

r GM 2l a 1 



M 2 = -[- 



;M 1; 



(4) 



vl l 2a 2 (1 + v 2 f/ 2 

where M\ is the rate of mass change of primary, M 2 is mass of secondary and a is the 
separation of two stars. The parameter a is taken to be 3/2 appropriate for Bondi-Hoyle 
accretion. The wind velocity of mass lost from primary is v w . This is difficult to determine, 
so it is simply taken to be the escape velocity from the surface, most likely an overestimate. 
The ratio between orbital velocity and wind velocity is v. 

Using the above, we calculate that a primodial 0.6 M Q star accretes 0.2 M & from the 
wind of a 7M star and the whole envelope of the primary is lost before thermal pulses 
begin. More detailed modelling and calculation of the secondary surface abundances are 
not carried out here. The abundances of the models are estimated for two extreme cases a) 
the accreted material mixes with the whole star and b) the accreted material remains in a 
surface layer. 

Th e abundances of C and N fit fairly well with the original data of IChristlieb et al 



(200 41). The abundances of C, N and O from the 3-D atmospheric analysis published by 



Collet. Asplund fc Trampedachl (120061 ) are significantly lower (by 0.8 dex or more) than 



the 1-D analysis (Table 7). Thus, it is possible that the observed high carbon and oxygen 
abundances of HE 0107-5240 are due to accretion from the wind of an AGB star. Compared 
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with oxygen abundances given by 



Bessell et al. 



(120041 ). oxygen is underproduced in the 



model. We investigate whether changes to the stellar ph ysics could result in a better fit to 
the oxygen. First we use the 12 C(a, 7) 16 reaction rate of lAngulo et al.l (119991 ) and find that 
increasing this rate by a factor of five is sufficient to raise the mass of 16 O dredged up to 
the obse rved value at neglig ible expense of the 12 C. However, the upper limit of the rate 
given by lAngulo et al.l (119991 ) is only 50 per cent greater than the adopted one in the relevant 



temperature range. Though some still doubt this result (llmbriani et al 



200 ll ) , there is as yet 



no evidence for an unknown resonance in the reaction, so the current reaction rate is too low 
to reproduce th e observation. Seco ndly we investigate the addition of convective overshooting 



as described in 



Pols et al. 



(119981 ) but find the effect small because of the stability of the 
Schwarzschild boundary at the base of the convective envelope during dredge up. Other 
possibilities include a higher oxygen abundance in low-metallicity stars compared to solar. 
This would imply the star has an initial composition enriched by the earlier generation of 
stars which is possible given [Fe/H] ~ —5 

Another notable difference is that our models have a surface enhancement of sodium 
by a factor of two which is due to the high intershell temperature of a 7M Q star. However 
sodium is observed to be enhanced by more than a factor of six. This can be fixed by either 
increasing the reaction rate of 22 Ne(p,7) which has an uncertainty of 10 or by lowering the 
reaction rate of 23 Na(p, a). Also, during thermal pulses, sodium is produced when 22 Ne from 
the intershell reacts with protons in the hydrogen shell. Sodium could then be created at the 
base of the convective envelope by hot-bottom burning and would be brought to the surface 
in a 7M Q star. The surface carbon to nitrogen ratio may also change through hot-bottom 
burning but it would still be within the observational range after few thermal pulses. The 
strength of thermal pulses is too weak for carbon or s-process elements to be dredged up. 
Alternatively, the over-abundances of sodium may again be due to enrichement by an earlier 
generation of stars. 

The second star, HE 1327-2326, has strong enrichment in the neutron-capture element 
strontium but no barium. Strontium is usually the product of the weak s-process, which 
produces essentially no barium. This would fit the binary scenario if the mass accretion 
continues briefly after the start of some third dredge-up, when the surface of the primary 



is enriched in s-process neutron-capture elements. However, the neutron flux has to 



je low , 



and the efficiency of the weak s-process at low metallicity is uncertain (lAoki et al. 



20061 ). 



Also our Z = models do not show third dredge up. Another possible explanation is that 
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Element 


1JJ log e^A. ) 


oU log e(A ) 


With mixing 


Without mixing 


c 


6.81 


5.72 


5.44 


6.08 


N 


4.93 


2.91 


3.10 


3.71 





5.66 


4.95 


3.49 


3.71 



Table 7. Observed abundances of HE 0107-5240 compared with binary models, log e(X) = logio(Nx/Nn) + 12.0 



HE 1327-2326 was formed f rom a cloud enriche d by a nearby supernovae which produced 



strontium via the r-process (ITruran et al 



20021 ) . HE 1327-2326 is also much more enriched 



in nitrogen which can be produced by hot-bottom burning in intermediate-mass AGB stars. 
On the whole, the abundances seem to suggest that binary mass transfer at a later AGB 
phase might fit the observed abundances of HE 0107-5240. Such stars, which require the 
solution of a number of complex numerical problems, will be modelled in future work. 



7 CONCLUSIONS 

Mass transfer to a companion of low mass during AGB phases results in chemically peculiar 
stars such as carbon-enhanced metal-poor stars. Mass accreted by the secondary can signif- 
icantly enhance the surface abundances of the CNO isotopes and possibly neutron-capture 
elements. For example, the abundances of HE 0107-5240 can be modelled by wind accre- 
tion from a 7M primary during early AGB phases. HE 1327-2326 could also be formed 
in similar way but with mass transfer occuring at a later stage. If such mass transfer is 
non-conservative metals are also released to the interstellar medium after second dredge up 
and enrich the interstellar medium. This is possible even if the binary is sufficiently close 
that Roche lobe overflow is followed by common envelope evolution in which the envelope is 
lost from the system. These binary systems contribute more to the interstellar medium than 
a single star at the same stage of evolution because the envelope is lost earlier. However, the 
total contribution from a primary star in binary system is smaller than in a single system 
when the single star ends its life as a type 1.5 supernova. This would release a large yield 
of carbon, oxygen, nickel and iron. The combined contribution of binary stars and single 
stars depends, not only on the binary fraction and separation in the early Universe, but also 
the mass-loss rate and evolutionary endpoints of single stars of zero metallicity. The more 
stars in close binaries, the smaller would be the total nitrogen yield because the hot bottom 
burning in the AGB phase would be curtailed. 
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